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Abstract 

We studied the Fermi and Gamow- Teller responses of cold symmetric nuclear matter 
within a unified dynamical model, suitable to account for both short- and long-range 
correlation effects. The formalism of correlated basis functions has been used to construct 
two-body effective interactions and one-body effective weak operators. The inclusion 
of the three-body cluster term allowed for incorporating in the effective interaction a 
realistic model of three- nucleon forces, namely the UIX potential. Moreover, the sizable 
unphysical dependence of the effective weak operator is removed once the three-body 
cluster term is taken into account. 



1. Introduction 

The understanding of the interaction of low-energy neutrinos with nuclear matter is 
required for a description of a number of properties of compact stars. 

The neutrino mean free path is thought to play a crucial role in the mechanism leading 
to supernovae explosion, while neutrino emission is the main process driving the early 
stages of neutron stars' cooling. 

The charged-current weak response of isospin symmetric nuclear matter has been an- 
alyzed within the framework of correlated basis functions (CBF) theory in Refs. [H, 0- 
Unlike the earlier applications of the CBF formalism to the density and electromagnetic 
weak response 0, U Q, the works of Refs. [H Q are based on effective interactions and 
effective transition operators, allowing for a consistent treatment of short-range corre- 
lations, which are known to be dominant at large momentum transfer, and long-range 
correlations, leading to the excitations of collective modes at low momentum transfer. 

The effective interaction of Ref. [H has been obtained at two-body cluster level, in 
the cluster expansion of the ground-state expectation value of an hamiltonian including a 
truncated version of the Argonne wis potential Q . The authors of Ref. also took into 
account the effects of interactions involving three or more nucleons through the density 
dependent phenomenological potential originally developed in Ref. 0. In this work 
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we improve the CBF effective interaction by explicitly including the three-body cluster 
contributions, which allows for a more realistic description of three-nucleon forces at 
microscopic level. In particular we have been able to include the leading contributions of 
the UIX three-body interaction . The effective transition operators arc also consistently 
evaluated including the leading order three-body cluster terms. 

An overview of the formalism is given in Section [51 In Section [3] the effective inter- 
action at three-body cluster level for an hamiltonian with Argonne Vq + UIX potential 
is derived. Section [J] is devoted to the development of the weak transition operator at 
three-body cluster level. The correlated Fermi gas (CFG), the correlated Hartree-Fock 
(CHF), and the correlated Tamm-Dancoff (CTD) approximations are discussed in Sec- 
tion [5l The numerical calculations of the response and of its sum rule are reported in 
Section [6l Finally, conclusions and future perspectives are discussed in Section [T) 



2. Formalism 

In the low-momentum transfer regime (|q| of the order of 10 MeV), the non relativistic 
limit of the weak-charged current matrix element is expected to be applicable. The 
nuclear response to weak probes delivering energy w and momentum q at leading order 
in |q|/m reads 

S{q, ^)^jYl l(*n|Oq|vI/o>|'<5(u; + Eo- E^) . (1) 

In the above equation A is the particle number and Oq the one-body weak operator 
that induces a transition from the ground state |\l/o) to the excited state l^'n), which are 
eigenstates of the nuclear hamiltonian with energies Eq and En , respectively, i.e. 

7?|*„) =S„|*„). (2) 

The non relativistic Fermi (F) and Gamow- Teller (GT) operators describing low en- 
ergy weak interactions are 

i i 

0^^ = Eor«=.9AE (4) 

i i 

where gv = 1-00 and gA = 1-26 are the form factors at zero momentum transfer, while 
is the isospin-raising operator acting on the i-th nuclcon. 

The spin-longitudinal and spin-transverse components of the Gamow- Teller response 
functions are defined as the components parallel and orthogonal to q, respectively. They 
can differ significantly at large values of |q| and, in principle, should be calculated sep- 
arately. However, whenever not otherwise specified, with "Gamow- Teller response" we 
will be referring to the total response, defined as the sum of the cartesian components 

5«^(q,c.)=i EK*"l^qIl*o>|''5(.. + i?o-i?„). (5) 
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The differences between the integrated spin-longitudinal and spin-transverse responses 
will be discussed in Section [Ql 

Within CBF, the states appearing in Eq. ([T|) are written in the form 

where |<f>„) is the Slater determinant of non interacting n— particle n— hole state. The 
structure of the correlation operator, reflects the complexity of the Argonne v'q 
nucleon-nucleon potential 111. Il2l. Il3| 

A 

= 5 n 4 ' (7) 

3>i=l 



ith 



4 = Er(^..)of, (8) 

p=i 

and 

In the above equation, aij = Si ■ Oj and Tij = Ti ■ Tj, where cfi and Ti arc Pauli matrices 
acting on the spin or isospin of the i-th, while 

^..-(3rSrf,-<5"Vr^f , (10) 

with a, /3 = 1, 2, 3, is the tensor operator. 

Note that the symmetrization operator S is needed to fulfill the requirement of anti- 
symmetrization of the state \'^n), since, in general, [Of^ , O^^.] ^ 0. 

The variational parameters determining the shape of the radial functions f^{rij) have 
been fixed in Ref . [IJ] , minimizing the variational ground-state energy. In this work we 
will use the results corresponding to the hamiltonians with Argonne v'q and Argonne v'q 
+ UIX potentials. 

Following Refs. [H, 0: we only consider transitions between the correlated ground- 
state and correlated Iparticle-lhole {Ip — Ih) excited states. The np — nh states with 
n > 2 give a smaller contribution, mainly at large excitation energy.the weak response. 

The CBF matrix clement between the ground-state and Ip — Ih excitation reads 

where Pm and hi denote the whole set of quantum numbers of the single nucleon state, 
namely the momentum, the spin and the isospin projections along the z— axis. 

This quantity, entering all our calculations of the response function, will allow us to 
define the effective weak operators, as discussed in Section 
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3. Effective interaction 



Using the formalism of CBF and the cluster expansion technique, the authors of Ref. 
[isl [l[ , were able to develop an effective interaction, obtained from the bare Argonne 
potential, which incorporates the effects of the short-range correlations. In Ref. [l6l |. 



the two-body effective interaction of Ref. [15|, [1| was improved with the inclusion of the 
purely phenomenological density dependent potential of Ref. 17|, accounting for the 



effects of interactions involving more than two nucleons. The CBF effective interaction, 
vl2'^ , suitable for use in Hartree-Fock calculations, is defined through the matrix elements 
of the hamiltonian in the correlated ground-state 



(^oK'il^l'fo) 



(12) 



As suggested by the above equation, the effective interaction allows one to calculate 
any nuclear matter observables using perturbation theory in the orthonormal FG basis. 
However, in general, extracting the effective interaction is a very challenging task, in- 
volving difficulties even more severe than those associated with the calculation of the 
expectation value of the hamiltonian in the correlated ground state. 

3.1. Two-body cluster 

The procedure developed in Ref. [ij consists in carrying out a cluster expansion of 
the Ihs of Eq. ((T^ and keeping only the two-body cluster contribution. The sum of the 
two-body cluster contribution of the potential and kinetic energies are readily found to 
be 



= '- / rfriaCTr 

2b 2 



F12V12FU (Vii^i2)(ViFi2) 

m 



X (1 - AT^L) 



(13) 



where the symbol "CTr" denoting the normalized trace over the spin-isospin degrees of 
freedom of particles 1, 2 and 3 The spin-isospin exchange operator and the Slater function 
are defined as 



ir^j) = 3 



s.ii\{kprij) — kpTij cosikprij) 
[kprijY 



(14) 



being kp the Fermi momentum of the system. 

On the other hand, the expectation value of the effective potential is given by 



(^0|«i2^^|<f>0) 



P 



drizCTr 



(15) 



Therefore, the effective potential at two-body cluster level turns out to be 



'hi 



= Fvi2F (Vi^^i2)(Vi^^i2). 

2b m 



(16) 
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The effective potential of the above equation shghtly differs from the one reported in 
the hterature. The authors of Refs. [H Q have not done the integration by parts leading 
to the kinetic term of Eq. (|13|) . As a consequence they have neglected the terms in which 
the gradient operates on both the correlation function and on the plane waves. In our 
effective potential, these terms, although small compared to the other contributions, are 
fully taken into account. 

3.2. Three-body cluster 

We have improved the effective potential by adding the tree-body cluster contributions 
to the energy per particle firstly calculated in Ref. [18[. This allowed us to consistently 
include in the effective interaction the UIX potential, whose leading order terms emerge 
at three-body cluster level. 

The three-body cluster contribution appearing in the expansion of J^twi2J^ is given 



Within the FR diagrammatic scheme, the three-body cluster contribution to (W12) it is 
not merely the expectation value of the latter result, unlike the two-body case. As a 
matter of fact, the reducible diagrams arising form four-body cluster term of J-^vi2J-, 
the detailed calculations of which can be found in Ref. [l9[, needs to be taken into 
account. 



by 




(17) 
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Figure 1: Four-body reducible diagrams, v[ 



dir 



and their three-body reduction. 
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The direct term of the three-body cluster contribution in the FR expansion scheme 
is given by 



dir p2 
3h ~~2 



driarfrisCTr 



(18) 



(5Fi2-Fl3-f23)i'l2('5Fl2^'l3^23) — -^^12^12^12 (^^13 + -^23 ~ 1) ■ 

It includes the term 

{S Fi2Fi^F23)vi2{S Fi2Fi^F2'i) — ^i2i'i2-P'i2 (19) 

from the three-body cluster contribution of F'^vi2F, whereas the reducible four-body 
diagrams of Fig. [1] contribute with the factor 



A2Wl2A2(A'3+^l3-2). 



(20) 




Figure 2: Four-body reducible diagram, v^^^^,},' ^^'^ three-body reduction. 



The sum of the diagrams where particles 1 and 2 arc exchanged gives 



"12 



= - ^ y dri2dri3^?2CTr 

{SFi2FuF23)vi2iSFi2Fi3F23)P^J - A2«12A2(A'3 + 4^3 " • (21) 



The corresponding four-body diagram producing the term |fi2Wi2-F'i2(A^3 +-^23 ^'^)Pi2 
is drawn in Fig. [5] 

The diagrams in which particles 1 and 3 arc exchanged contributes with 



^^12 



^ / dri2dri3^^3CTr 



(5Fi2Fi3F23){)i2(5Fi2Fi3F23)P{^3^ — Fi2'0l2^12-F'f3-F'r^ 



S-' 13 



(22) 



where the term \Fi2Vi2Fi2{F'l^ - 1)A'3^ 

comes from the four-body reducible diagram of 

Fig. m 
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Since the potential is invariant under Xi X2- the diagrams with the exchange 
between particles 2 and 3 give the same contribution reported in Eq. ((22|) . The associated 
four-body reducible diagram is very similar to the one of Fig. [3]but with the loop attached 
to particle 2 instead of particle 1. 

Consider the diagrams with the circular exchange involving particles 1, 2 and 3. In 
this case there are no reducible four-body diagrams that partly cancel the reducible part 
of the three body diagram. In addition, there are no three-body reducible diagrams with 
circular exchange at all. However, the four-body diagram of Fig. [31 with no correlation 
lines linking particles 1 and 2 to the others, can be reduced to a three-body term, so that 
the three-body diagram with a circular exchange reads 



3h 



dri2dri3^i2^i3^23CTr 



[S Fi2Fi2,F2'i)vi2{S Fx2FizF2z)Pi2 Pvi " Fi2Vl2Fi2Ff.^P1^ PI 



(23) 



The three-body cluster contribution to the PB kinetic energy contains terms of the 
kind y\{S Fi2FizF2z) ■ Their explicit expressions can be obtained from the corresponding 




-1 X 




Figure 3: Four-body reducible diagram, v^?-^^,, and its three-body reduction. 
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Figure 4: Four-body diagram, f|J,'l_>3i,, 
exchange between particles 1, 2 and 3. 



that contributes to the three-body diagrams having a circular 
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equations for the two-body potential by substituting the first term of the normalized 
traces with 1181 



Vl2iSFi2Fi3F23) ^ - 2[5(V?Fi2)A3^^23] - 2[5(ViFi2) • (Vi^^i3)F23] , (24) 



while 



V12F12 ^ -2(V2Fi2) 



(25) 



for the second term. Terms with (VfFi2) are denoted by W^™ , those having (V1F12) • 
(V1-F13) are included in U. On the other hand, the three-body cluster terms belonging 
to Wp arise from the diagrams where particles 1 and 2 are exchanged 



(f) 



Wf 



dri2dri3£i2^'i2^i2 • CTr 



(5A2A3A3)[5(ViFi2)A3A3)A7 - A2(ViA2)(A'3 + - 



(26) 



and from the ones with circular exchange 



(T) 



ctr 
3b 



-P 



dri2dri3£i3^23^'i2n2 ' CTr 



(5A2A3A3)[5(ViA2)A3A3]A7A7 - A2(ViA2)A'3A7A7 



(27) 



The contributions to U stem from the diagrams with the exchange P13 



(T) 



Uf 



^—p' / dri2dri3f (ri3)fi3 • CTr 
36 m 



( 5Fi 2 i^i 3 F2 3 ) [5 ( V 1 i^i 2 ) Fi 3 F2 3 ) -P] 
and from those having circular exchange 



(28) 



(T) 



Uf 



ctr 
35 



dri2dr 



13*12*13 



^'13^23 • CTr 



{SFi2Fi3F23)[S{WiFi2)Fi3F23]P^2"Pl 



(29) 



Note that in this case there are no subtraction terms arising from reducible diagrams. 

Finally, the three-body cluster term of the three-body potential V123, like UIX model, 
can be easily obtained by keeping only the first term in the traces of Eqs. (fT5)) . (PT|) . 
(|22|) and ((23l) with the following replacement 



V12 -> 



123 



(30) 



As in the construction of the density dependent potential from UIX, developed in 
Ref. the issue of the exchange pattern has to be carefully analyzed. The distinctive 
feature of the present calculation is that vl2^ contains the correlation between particles 
1 and 2 making possible to implement the inversion of P^"^ in a straightforward way. To 



be definite, consider the three-body cluster contribution of the ground-state expectation 
value of the two-body potential 



^'12 



3b 2 



(31) 
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Figure 5: Central (a) and the tr (b) channels of the effective potentials at two-body and three-body 
cluster level calculated for SNM at p = 0.16 fm^^. 



Note that the above equation summarizes the terms corresponding to Eqs. (|18p . (HU, 
and ([^5]) . A comparison with Eq. (fTS]) immediately leads to 



-eff 



3b 



=p / drgCTrs ^(xi, .T2; 2:3) 1 - Fj^^fa 



2Pi7^^3 + 2Pr/Pi7^i2^i3^23 ) (1 - Pi7el2r' 



(32) 



where the subscript 3 indicates that the spin-isospin normalized trace has to be performed 
over coordinates of particle 3 only. A similar argument holds for the three-body cluster 
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Figure 6: EoS of SNM in the low density regime. In the panel (a) the bare hamiltonian only contains the 
Argonne v'g potential, while in the panel (b) the UIX three-body interaction model is added to it. The 
dotted and the solid lines display the three-body and the two-body cluster effective potential results, 
respectively. The FHNC/SOC calculations are represented by the shaded region, accounting for the PB 
and JF kinetic energy difference. 



term of the kinetic energy and the three-body potential contributions to the effective 
potentiaL 

In Fig. [S] the central and tr components of the effective potentials at two-body and 
three-body cluster level are compared. The starting bare NN interaction is the Argonne 
v'q; for the three-body cluster results the UIX three-body potential has been included in 
the calculations. 

Starting from a bare hamiltonian only containing the Argonne v'g NN potential, we 
have computed the EoS of SNM for the low-density regime using both the new three-body 
cluster effective interaction and the older one with only two-body cluster diagrams. The 
results have been compared with the corresponding FHNC/SOC calculations, displayed 
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as a shaded region Fig. [6] to account for the PB and JF kinetic energy difference. The 
curve corresponding to 'y^2''^|gj^ is much closer to the FHNC/SOC results than the one 

obtained with the older fi2'^L • 

In the lower panel of Fig. [61 the EoS of SNM are shown for an hamiltonian containing 
the Argonne Vq NN potential along with the UIX three-body interaction model. Again 
the curve obtained from the three-body cluster effective potential is close to the full 
calculation and, it exhibits the saturation, which is a remarkable feature. 

As a second step in the definition of the effective interaction, we have adjusted the 
variational parameters of the correlation functions in order for the effective hamiltonian 
exactly reproduces the energy per particle obtained with the full FHNC/SOC calculations 
at saturation density. In particular the "healing distance", of both the central, dc, and 
the tensorial, dt , correlations and the quenching parameters ap [3l , which can be given 
the interpretation of the low-energy parameters of the effective interaction, have been 
reduced. 



4. Effective weak operators 

In the case of Ip — Ih excitation, the effective weak operators O^^^ , introduced in 
Ref. [H, are defined through the relation 

For the sake of giving a unified description of matrix elements associated with both 
the Fermi and Gamow- Teller transitions, it is convenient to distinguish the common 
radial part from the specific zero momentum form factors and spin-isospin operators, 
defining 

Oq(l) = e*'i--id,,(l) . (34) 

From Eq. (g]) it follows that 6l^{\) = gyr^ and 6^7(1) = QA^irt ■ 

As for the calculation of the hamiltonian expectation value, a cluster expansion of 
the weak operator correlated matrix element can be performed [lH |. The smallness 
parameters in this case are f![j — 1 and ffy Following Ref.Q, we denote b y cj , the 
quantum numbers of occupied states in both $o a-nd It can be shown [3l.ll9| that 

the CBF matrix element of the effective weak transition operator takes the form 



The term X^cC**) contains all the connected diagrams with one hi vertex or one hi- 
exchange line. An analogous definition applies for ^(j{hi). On the other hand, ^(^{q',Prrn hi) 
amounts of connected diagrams having one single Pm^ii-exchange line or one single Pmht 
vertex. Moreover, the weak operator Oq(l), which carries a momentum q and a spin- 
isospin operator, is attached to the point 1 of the diagrams belonging to '^c(ci]Prm hi). 
Note that only c,; states are present in both the bare vertex and the bare exchange lines. 
In other words, unlike in the cluster expansion for the energy per particle, the hole state 
hi is lacking. 
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Figure 7: Non vanishing three body diagrams at first order in / — 1 emerging from the numerator of Eq. 
ll35ll . 



The two-body cluster diagrams coming from both the numerator and denominator of 
Eq. (|35p have been computed by the authors of Ref. [l^ , while in Ref . 0, [23| a different 
truncation scheme has been adopted and the numerator has been approximated to the 
unity. 

In this work, in addition to the numerator and denominator two-body cluster di- 
agrams, for the sake of consistency with the effective interaction, we have calculated 
the three-body cluster diagrams at first order in / — 1. They are associated with the 
numerator of Eq. ([55]) , as they originate from 

{/23-l,Oq(l)} =26q(l)(/23-l). (36) 

The only non vanishing three-body cluster diagrams are depicted in Fig. [T] In the 
thermodynamic limit, diagram N3bia reads 

X ^(ai«2ap„|a.(l)(/23 - l)A7l«ia2«/..) , (37) 

where \ai) denotes the spin-isospin state of particle i. The discretized momentum con- 
servation is expressed by the Kronecker delta function. By computing the spin-isospin 
matrix elements, whose values are reported in Appendix E of Ref. |19j . one can show 
that in the limit of zero momentum transfer, q — > 0, A^3&ia cancels the contribution of 
diagram N2hia, represented in Fig. [8] (denoted as Fid j and GTld j in Ref. [Hj). This 
is an indication that three-body diagrams need to be taken into account and they play 
a relevant role in the sum rules of the weak response, which will be estimated at a later 
stage. 

It can be readily shown that the analytic expression of diagram NSbib is given by 

X J2{aia2apjd^ril){f23 - l)Pr2"A7|aia2«h.) • (38) 



The evaluation of the spin-isospin matrix elements can again be found in Appendix E of 
Ref. [H. 

12 



Pmhi 

. 
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Figure 8: Two-body diagram of the first order term in / — 1 coming from the numerator of Eq. II35II . 

5. Response functions in different approximations 

5.1. Correlated Fermi gas (CFG) and correlated Hartree-Fock (CHF) 

In both the correlated Fermi gas (CFG) and correlated Hartree-Fock (CHF) approx- 
imations, the weak response of cold SNM, defined in Eq. ([Ij, is given by [l[ 

Spcici.u:) -1 \i%..-M\Ot!^m?5{'^ + ^p,^~^K)- (39) 

Vm hi 

Within the CFG approximation, the single particle energies are those of the non-interacting 
hamiltonian 

e„. = |^. (40) 

The single particle approximation is retained in the CHF; however the potential enter 
the calculation of the single particle energies. It has long known [2lj that the Hartree- 
Fock approximation is not suitable for nuclear potentials having a repulsive core, like the 
Argonne models, because it does not encompass the correlations between nuclcons. We 
use instead the effective potential described in Section [3l which is appropriate for mean 
field calculations. The single particle energy is then given by 



1^2 A 



(41) 



where A is the antisymmetrizing operator. The single particle wave functions are plane 
waves 



where rja^ = XciXri represents the product of Pauli spinors describing the spin and the 
isospin of particle i and V is the normalization volume. 

While in CFG calculations the correlations enter only through the effective weak op- 
erators, within the CHF approximations they sizably determine the effective interaction. 

In the case of SNM {v = 4) for potentials of the form of Argonne wis, carrying out 
the summation over the occupied states with |kj| < kp yields 

+ / ["^^- ~ + + + ■ (43) 

For the sake of simplicity, in the latter equation the superscript "e//" has been omitted. 
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5.2. Correlated Tamm-Dancoff (CTD) 

Since the correlated Ip—lh states are not eigenstates of the fuh nuclear hamiltonian, 
transitions between them are in principle allowed. They can be accounted for within the 
Tamm-Dancoff approximation, that amounts to expanding the final state of Eq. ^ in 
series oi Ip — Ih excitations. 

Because the hamiltonian is translationally invariant, the total momentum q of the 
state is conserved, and the momenta of the particle, Pm, and the hole, h^, satisfy the 
relation Pm — = q. 

The nuclear hamiltonian commutes with the total isospin, T with the total isospin 
projection along the z— axis, T^, and with the total spin, S. However, because of the 
tensor term of the potential, the hamiltonian does not commute with Sz, the total spin 
projection along the z— axis. 

The combinations of particle hole pairs that are eigenstates of S and Sz , that define 
the particle-hole Clebsch-Gordan coefficients, are shown in Table \T\ The differences 
between the total spin states of the particle particle pairs, also given in Table [H are due 
to the phase factor appearing in the canonical transformations to particles and holes [2^ . 
The treatment of the total isospin can be done in complete analogy, replacing the up 
and the down single particle spin states with the proton and the neutron isospin states, 
respectively. 

Table 1: Spin configurations for a particle particle pair and a particle hole pair for spin- 1/2 particles. 



Total spin state particle particle particle hole 



5* = 


l,Sz 


= 1 


tt 




-n 


5* = 


l,Sz 


= 




-;t) 






l,Sz 


= -1 








S = 


0,5, 


= 




-;t) 





It is possible to reduce the computational effort necessary for solving the Tamm- 
Dancoff equations, by starting with combinations of with definite T, Tz, S and 

\^n)^?.i- E C;:^:'''\%,..^.)tt.ss., (44) 

A further simplification arises by noting that the final states of both the Fermi and the 
Gamow- Teller transitions are characterized by having T ~ 1 and = 1 . To simplify the 
notation, the isospin indexes may then be omitted 

l<i>n>?^^= E Cpfh;l'i>P™;h.>55., (45) 
and it is understood that T = 1 and = 1. 
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The eigenvalue equation for the effective hamiltonian defines the excitation energy 



H^"\^n)r^ = ( E + E ^f) 1*")^^ = (Eo + )|$„>?^^ . (46) 

i i<j 

Multiplying from the left the previous equation by 55^ (*i'p„:h | and using the orthonor- 
niality of the Ip — Ih states yields 

E 55.(<i>p„;h,|^'^^^|<i>p„,;h,>55^ = C;f^; (Eo + U^^) ■ (47) 

Thus, finding the coefficient Cp amounts in diagonalizing the block diagonal hamil- 
tonian for the two subsets of the Ip ~ Ih basis having T = 1, Tz = 1 corresponding to 
S ^ and to 5* = 1 . This is a much less expensive computational task than diagonaliz- 
ing the hamiltonian in the full Ip — Ih basis. As a consequence, this approach allows for 
considering a larger number of momentum states. 

It has been shown [l^ that, singling out particles 1 and 2 from both the ground state 
and the Ip — Ih excited state, the matrix element of the hamiltonian reads 

ss,($p„;hji?"'^^|$p,„;h.)ss_i = (^o + ep,„ - eh,)Sp,,^pJh,h^5s,s', 

+ {Pnh,\vli^\hjp^)ss'^s.] (48) 

As far as the notation is concerned, with {pnhj\0i2\hiPm) we denote the two-body 
matrix element of the operator O12 



(Pn hi\Oi2\hj Pm) = J dxi,2(l}l^{xi)4>l.{x2)Oi2A[4'hj{Xl)(l)p,„iX2)]- (49) 

In the two-body matrix element of the effective potential (p„ hi|{i^2''^|hj Pm)ss',s^ the 
spin and isospin projections along the z axis of the particle hole pairs Pm — hi and p„ — hj 
are combined as in Tablc[T]to have definite S, and total spin projections along z equal to 
S'z and 5*2, respectively. We recall that the total isospin and its z-projection are T = 1 
and Tz = 1. 

The direct and exchange terms for the case 5 = 0, relevant to the Fermi transition, 
for SNM read 

(p„h,|{)^|-^|h,p,„)[joo = -^ / rfri2e-"i--- i;[2 

(p„ h,|i;i2|h, p„)§oo = ^ / dri2e'''--^^nz;j2 " ^^[2 + " 3<J) , (50) 

where k^j = — hj. Again, for the sake of simplicity, the superscript ^^eff" has been 
omitted where the channels of the effective potential are specified. 

For the Gamow Teller transition, the final state has S* = 1; hence it is necessary to 
compute the nine matrix elements (p„ 15^2''^ |hj p,„) 55/ 5, corresponding to Sz,S'z = 
-1,0,1 
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(p„ h,\vll^ \hj Prr,)^^:^ 
(p„ h,\vl{,^ \hj p,„)l^^ 

(p„ h,\vl(^\hjp,n)'ioi 
(p„ h,|?)i^^|hjP,„)ioi 
(p„ hi\vl^^\hjp„^)i_ll 
(p„h,|z)i^-''|hjp,„)'j_ii 
(p„ h,\vf/\hjp„,)%Q 
{Pn hi|w;=^^|hj P™)^oo 



2 
1 



2<J + z;- 1 



3^12 



+ K2--^2)(l-3^ 



12 

J ^12 



J ' 1 



(a;i2 -jyi2f 

2 



dri2e-*i'"" 



-"12 "^12 



'12 



12 



= - / dri2e*''---- 



2K2-«r2)(l-3^ 

^ '12 



^^12-^^12-^12+^^12 
2 



(51) 



Replacing the final state of Eq. ([^S]) in the definition of the response, Eq. ([T]) yields 



(52) 



6. Numerical calculation of the response 

We model the infinite system using a cubic box of side L = V^^^ with periodic 
boundary conditions. Hence, the single particle wave functions are the plane waves of 
Eq. ((42l) with the discrete momenta 



nt, =0,±1,±2,. 



(53) 



For zero temperature SNM, all single-particle states with |k| < kp are occupied in the 
ground state. The momenta of the \p — Ih excitations are such that |hj| < kp and 
\Pm = hi + q| > fci?. For the hole and particle momentum to be on the lattice of allowed 
momentum states in the box, the momentum transfer must be such that 



2n 



Ug^ = 0,±1,±2, 



(54) 
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Given the magnitude and the direction of the momentum transfer, the side of the box 
can be determined by inverting the latter equation 



Ztt I 

^=|^V<+<+<- (55) 



The size of the basis, that can be increased by increasing ^Jn^^ + n^^ + n^^ , has been 

determined requiring that the response of a system of noninter acting nucleons computed 
on the lattice agreed with the analytical result of the FG model P, 0] ■ The FG response 
is obtained replacing the effective operator with the bare operator in Eq. p9p and using 
the single particle energy of Eq. ([^01) . 
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Figure 9: Fermi response functions calculated at q = 0.3 fm ^ . Panel (a) shows the response calculated 
at two-body cluster level for both Oq^^ and f ^2^ different choices of correlation functions. In panel 
(b), the three body cluster has been included. The plus marks show the response for a non interacting 
FG. 

The analytical calculations can be performed using a continuum of momentum states 



22], while the numerical result consists of collection of discrete delta function peaked 
at the values of the single particle energies. For a better representation of the results, 
as well as for fitting purposes, a gaussian representation of the energy conserving delta 
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function has been adopted 



S{x- .To) 



cxp 



(56) 



For sufficiently small values of the gaussian width cr, the results become insensitive to 
it. All the results that will be shown in this section refer to SNM at equilibrium density 
p = 0.16 fm~^. 




Figure 10: Same as Fig. |9]but for Gamow- Teller response functions. 



6.1. CFG 

FHNC/SOC calculations and their associated minimization procedure provide a set 
of correlations function, corresponding to the minimum of the liamiltonian expectation 
value. We have found the best correlation functions for the Argonne Wg, Wg two-body 
potentials, and for comparison we have also considered the correlations of Ref. [23| 
corresponding to Argonne wig. With these correlations, the Fermi and Gamow- Teller 
response functions have been evaluated in CFG approximations with the sake of testing 
the contribution of the three-body cluster in the effective weak transition operators. The 
effective potential indeed does not enter the CFG calculations. 
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When only two-body cluster diagrams are considered, as in Refs. P, il, the CFG 
response, suppressed by a 20 — 25% with respect to the FG case, exhibits a sizable 
dependence on the choice of correlations, as shown in the panels (a) of Figs. [5] and [TU[ 
These figures refer to a transfer momentum q ~ |q|(4a; + Ay + 4z)/-\/48 with |q| = 0.3 
fm"-'^. The folding gaussian function has a width of 0.25 MeV. 

This unphysical effect is removed once the effective weak transition operator is com- 
puted at three-body cluster level. As a matter of fact, the CFG curves in the panels 
(b) of the aforementioned figures are very close, when not superimposed, to each other. 
Therefore our results appear to be more robust than those of Refs. P, [lli as the phys- 
ical quantities should not be sensitive to the details of the short range behavior of the 
correlation functions. 



6.2. CHF and CTD 
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Figure 11: Fermi response functions calculated at g = 0.3 fm~^ using v^^i and iig/ + UIX interaction 
models. Panel (a) shows the full response across all the values of u); panel (b) is a magnification of the 
small u) region. Responses are folded with a Gaussian of width 0.25 MeV. 



The nuclear matter response calculated in CTD and CHF approximations for |q| = 
0.30 MeV is displayed in Figs. [TT] and [T^ for Fermi and Gamow- Teller transitions, 
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respectively. By comparing with the findings of Refs. [l|, it can be noticed that the 
peak corresponding to the collective mode is shifted to lower energies when the tree-body 
cluster is included. This effect, due to the change of the single particle energies, is only 
slightly mitigated when the UIX potential is included in the hamiltonian. 




Figure 12: The same as Fig. Illl but for Gamow- Teller response functions. 

Moreover, the three-body cluster produces a small depletion of the Fermi resonance 
at |q| = 0.30, more apparent when the nuclear hamiltonian is lacking of the three-body 
potential. 

It is worth remarking that the two-body cluster results of Ref. [l[ have been obtained 
using the first six operators of the Argonne potential rather than the Argonnc v'q used 
in the present work. However, we have performed two-body cluster calculations with 
Argonne v'q and the results are in extremely close to those of the truncated version of Wg . 

In Fig. [131 where the Fermi and Gamow- Teller responses are plotted for different 
values of |q| ranging from 0.10 fm~^ to 0.50 fm~^. The position of the highest peak of 
the Fermi responses shifts from |q| ~ 0.40 fm~^ of the two-body cluster approximation 
to |q| ~ 0.30 fm~^ of the present calculation. For the Gamow- Teller case, the three-body 
cluster contribution does not produce an appreciable shift, being the highest resonance 
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remains peaked around |q| ~ 0.25 fm ^. 



0.08 




Figure 13: Fermi (upper panel) Gamow- Teller (lower panel) response functions calculated at q 
0.10, 0.15, 0.20, 0.25, 0.30, 0.40, 0.50 fm"! using v^i+mX potential and correlations. 



6.3. Sum rules 

The set of final states in Eq. ([T]) is not exhausted by Ip— Ih exeitations. In principle, 
transitions to more eomplex multi p ~ h states should also be considered. So far, the 
contribution of these states have been neglected, however, an estimate of their importance 
can be obtained computing the sum rules. The static structure function is defined by 

5(q) = j dujSi(i,uj) 



-(*o|6J,Oq|vI'o) 



(57) 



21 



While a direct integration of the CTD response functions allows for the evaluation of 
S'(q), from the last line of the latter equation it is clear that the static structure function 
can be also evaluated by computing the ground state expectation value of OqOq. Within 
the correlated basis function approach, the variational ground-state (VGS) expectation 
can be expressed in terms of the two-body operatorial distribution functions, calculated 
by means of the FHNC/SOC scheme in Ref. [H. 

While the VGS calculations include all the multi p — h excitations in the CBF basis, 
only the correlated Ip ~ \h states are taken into account in the CTD approximation. 
Therefore, multi p — h contributions can be estimated from the difference S^'~'^{q) — 

Note that an interplay between many-body correlations and multi p ~ h excitation 
could in principle take place. As a matter of fact, while VGS includes many-body cor- 
relations through the chain summations, in the CTDA of Ref. [l| only two-body cluster 
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Figure 14: Static response function for Fermi, panel (a), and Gamow- Teller, panel (b) transitions. 
Results are shown for a noninteracting FG (solid lines), for the integral of the CTD approximation of 
^(q, Lo) (crosses) for the full FHNC/SOC VGS calculations (dashed lines) and for the three-body cluster 
approximation of the VGS (dotted-dashed lines). 
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terms have been considered. In the present work, three body cluster contributions have 
been taken into account in CTD calculations. Moreover, for the sake of isolating the 
effect of the multi p — h excitations from the one of the many-body cluster contribu- 
tions, we have also calculated the two-body operatorial distribution functions entering 
the variational ground state estimate of S{q) at three-body cluster level. 

The static structure functions for the Fermi and Gamow- Teller transitions are dis- 
played in panel (a) and (b) of Fig. [Ml respectively. Besides the FG curve, all the 
results have been obtained with an effective interaction that incorporates the the Ar- 
gonne Wg-f-UIX potential and the corresponding correlations are used in the calculation 
of the effective weak-transition operators. 

The curves corresponding to the Fermi transition are normalized in order for the 
sum rule of the non interacting FG to approach unity for large values of the momentum 
transfer. On the other hand, the Gamow- Teller results are normalized in such a way that 
both the transverse and longitudinal sum rules, to be defined in the following, tend to 
one for large |q|. 

As noted in Ref. the VGS results obtained by computing the two-body distri- 
bution functions with FHNC/SOC summation scheme, because of the approximations 
involved in the calculation, do not fulfill the condition S'(O) = 0, required by baryon 
number conservation. The three-body cluster variational results, denoted by VGSsf,, 
also does not fulfill the baryon number conservation; however the reason for this most 
probably lies in the three-body cluster approximation which is known not to fulfill the 
sum rules of the two-body distribution functions. Conversely, the static structure func- 
tion obtained within CTD approximation does exhibit the appropriate low-momentum 
limit. 

As far as the multi p — h excitations are concerned, the two-body cluster results of 
Ref. [l| show that their contribution is smaller than the dominant Ip — l/i excitation, 
but it is not negligible. When three-body cluster is accounted for, the VGS and the 
CTD curves get closer. The shift turns out to be detectable for the Fermi transition 
case, while for the total Gamow- Teller response is very small. Considering also the fact 
that VGSsb results lie well above the CTD crosses, we may conclude that the difference 
between variational ground state and correlated Tamm-Dancoff results has largely to be 
ascribed to the multi p — h excitations. 

There are experimental and theoretical indications that the spin longitudinal and 
spin transverse response functions can differ significantly due to tensor forces. Thus, 
we studied how the UIX three-body force affects these quantities computing the spin 
transverse and spin longitudinal static structure functions, defined as 

•' n 

S^q) ^jf dujJ2mM ■ d^^\^o)\'S{cu + Eo- E^) (58) 

A comparison between Fig. [m and Fig. 17 of Ref. d shows that the inclusion of the UIX 
potential brings the CTD curves for S'^{q) closer to those of VGS across all the values 
of |q|. For relatively large momentum transfer it is clear that the difference between the 
variational ground-state results and the CTD calculations is due to 2p — 2h excitations. 
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The curves VGS and VGSsf, are almost superimposed in this region, indicating that the 
contribution of many-body clusters is very small. 

As far as the longitudinal static response function is concerned, the position of the 
maximum of the CTD calculations including UIX potential almost coincide with the 
VGS and VGSab results. At small momentum transfer however, the three-body cluster 
CTD points lie below the variational results, mostly because of deficiencies of the latter, 
unable to fulfill the sum rules. 




Figure 15: Spin longitudinal (left panel) and spin transverse (right panel) static response functions, cal- 
culated using VGS (solid lines), two-body CTDA (squares), two-body CTDA (stars), and noninteracting 
FG (dotted lines). 

Recently [25l | , the sum rules relations (additional sum rules with an increasing power 
of u! in the integrand can be defined) has been inverted to obtain the spin-response 
function of PNM at zero momentum transfer. The authors of Ref. [1^ used the formalism 
of AFDMC to compute the operatorial distribution functions, obtaining very promising 
results. As a follow up of the present work, we are planning to compute the response 
function of PNM and compare it with their findings. 



24 



7. Conclusions 

Wc have shown that the weak response of isospin symmetric nuclear matter, obtained 
from the effective operators computed in the two-body cluster approximation, exhibits 
a sizable dependence on the choice of the correlation function, which is in fact unphysi- 
cal. When three-body clusters diagrams are considered, the transition matrix elements, 
in which correlations enter only through the effective weak operators, become nearly 
independent of the correlation functions. 

The three-body cluster contributions, described in Section [31 have been consistently 
included in the construction of the effective interaction. As a first step we have computed 
the EoS of SNM for an hamiltonian including the two-body potential only. In this case, 
the three-body cluster effective interaction provides a EoS of SNM much closer to to the 
one resulting from full FHNC/SOC calculations, cornpared to the one obtained using 
the two-body cluster effective interaction of Ref. [l|, [l^l . The leading contributions of 
the UIX potential, emerging at three-body cluster level, have been also included in the 
effective potential. As a result the EoS of SNM exhibits saturation at p ~ 0.18 fm^^. 

Inclusion of the three-nucleon interactions also affects the single particle spectrum, 
leading in turn to a shift of the CHF response as a function of energy transfer. 

The main effect of the three-nucleon force on the TDA response originates from a 
change of the off-diagonal elements of the effective interaction. As a result, the collective 
mode associated with the Fermi transition at |q| = 0.3 MeV turns out to be shifted 
to lower energy, although its magnitude is nearly unaffected by the three-body cluster 
contributions. On the other hand, a depletion of the peak is observed for the Gamow- 
Teller transition. The analysis of the TDA response for different momentum transfer also 
reveals a sizable effect of three-nucleon cluster contributions. 

The sum rules for the Fermi transition comes closer to the variational results once 
the three-body cluster cluster is taken into account, thus confirming the importance of 
many body effects, which are included in the variational calculations through the chain 
summations. The residual discrepancy cannot be accounted for by the n > 4— body 
cluster contributions, and is likely to be ascribable to the effect of multi p — h excitations, 
which are taken into account in variational calculations. This effect appears to be even 
larger in the structure function obtained from the Gamow- Teller response. In this case 
the change due to the three-body cluster is in fact very small. 

Some improvements are observed in the sum rules of the longitudinal and the trans- 
verse response. As shown in Fig. [151 the results of the three-body cluster calculations 
of S^{q) carried out within TDA are slightly closer to the variational ones for all the 
values of |q|, compared to the two-body cluster case. Moreover, the position of the max- 
imum of the TDA calculations of the static transverse response is almost coincident with 
that obtained from variational calculations. At small momentum transfer however, the 
three-body cluster contributions move the TDA S'^{q) away from the variational result. 

References 

[1] S. CowcU, V. R. Pandharipandc, Phys. Rev. C 70 (2004) 035801. 

[2] O. Benhar, N. Farina, Physics Letters B 680 (2009) 305 - 309. 

[3] S. Fantoni, V. Pandharipande, Nuclear Physics A 473 (1987) 234 - 266. 

[4] A. Fabrocini, S. Fantoni, Nuclear Physics A 503 (1989) 375 - 403. 

[5] A. Fabrocini, Phys. Rev. C 55 (1997) 338-348. 

25 



[6] R. B. Wiringa, V. G. J. Stoks, R. Schiavilla, Phys. Rev. C 51 (1995) 38-51. 
[7] I. Lagaris, V. Pandharipande, Nuclear Physics A 359 (1981) 349 - 364. 

[8] B. S. Pudliner, V. R. Pandharipande, J. Carlson, R. B. Wiringa, Phys. Rev. Lett. 74 (1995) 4396- 
4399. 

[9] J. W. Clark, P. Westhaus, Phys. Rev. 141 (1966) 833-857. 
[10] S. Fantoni, A. Fabrocini, in; J. Navarro, A. Polls (Eds.), Microscopic Quantum Many-Body Theories 

and Their Applications, volume 510 of Lecture Notes in Physics, Springer Berlin / Heidelberg, 1998, 

pp. 119-186. 10.1007/BFb0104526. 
[11] V. Pandharipande, Nuclear Physics A 174 (1971) 641 - 656. 
[12] M. L. Ristig, W. J. T. Louw, J. W. Clark, Phys. Rev. C 3 (1971) 1504-1513. 
[13] V. Pandharipande, Nuclear Physics A 181 (1972) 33 - 48. 

[14] A. Lovato, O. Benhar, S. Fantoni, A. Y. lUarionov, K. E. Schmidt, Phys. Rev. C 83 (2011) 054003. 

[15] S. Cowell, V. R. Pandharipande, Phys. Rev. C 67 (2003) 035504. 

[16] O. Benhar, M. Valli, Phys. Rev. Lett. 99 (2007) 232501. 

[17] I. Lagaris, V. Pandharipande, Nuclear Physics A 334 (1980) 217 - 228. 

[18] J. Morales, V. R. Pandharipande, D. G. Ravenhall, Phys. Rev. C 66 (2002) 054308. 

[19] A. Lovato, Ab initio calculations on nuclear matter properties including the effects of three-nuclcons 

interaction, Ph.D. thesis, SISSA-ISAS Trieste, 2012. 
[20] N. Farina, Weak Response of Nuclear Matter, Ph.D. thesis, Sapienza Universita di Roma, 2009. 

arXiv 0901.2507 
[21] R. Jastrow, Phys. Rev. 98 (1955) 1479-1484. 

[22] A. Fetter, J. Walccka, Quantum Theory of Many-Particle Systems, Dover Books on Physics, Dover 
PubUcations, 2003. 

[23] A. Akmal, Variational studies of nucleon matter with realistic potentials, Ph.D. thesis. University 

of Illinois at Urbana-Champaign, 1998. 
[24] V. R. Pandharipande, J. Carlson, S. C. Piepcr, R. B. Wiringa, R. Schiavilla, Phys. Rev. C 49 (1994) 

789-801. 

[25] G. Shcn, S. Gandolfi, S. Reddy, J. Carlson (2012). arXiv llSOS ■ 5^ 



26 



